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WAS STAR-FORMATION SUPPRESSED IN HIGH-REDSHIFT MINIHALOS? 
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ABSTRACT 

The primordial gas in the earliest dark matter halos, collapsing at redshifts z ^ 20, with masses Mhaio '--^ 10^ M©, 
and virial temperatures Tvir < lO'^K, relied on the presence of molecules for cooling. Several theoretical studies 
have suggested that gas contraction and star-formation in these minihalos was suppressed by radiative, chemical, 
thermal, and dynamical feedback processes. The recent measurement by the Wilkinson Microwave Anisotropy 
Probe (WMAP) of the optical depth to electron scattering, t « 0.09 ± 0.03, provides the first empirical evidence 
for this suppression. The new WMAP result is consistent with vanilla models of reionization, in which ionizing 
sources populate cold dark matter (CDM) halos down to a virial temperature of = lO^K. On the other hand, we 
show that in order to avoid overproducing the optical depth, the efficiency for the production of ionizing photons 
in minihalos must have been about an order of magnitude lower than expected and lower than the efficiency in 
large halos that can cool via atomic hydrogen (Tvii- > lO^K). This conclusion is insensitive to assumptions about 
the efficiency of ionizing photon production in the large halos, as long as reionization ends by z = 6, as required 
by the spectra of bright quasars at z <^ 6. Our conclusion is strengthened if the clumping of the ionized gas evolves 
with redshift, as suggested by semi-analytical predictions and three-dimensional numerical simulations. 



1. INTRODUCTION 

How and when the intergalactic medium (IGM) was reion- 
ized is one of the long outstanding questions in astrophysical 
cosmology, holding clues about the onset of structure formation 
in cold dark matter (CDM) cosmologies, and the nature of the 
first generation of light sources (see Barkana & Loeb 2001 for 
a review). The reionization history can be probed via the scat- 
tering of the cosmic microwave background (CMB) photons 
by the free electrons liberated during reionization (see, e.g., 
Haiman & Knox 1999 for a review). The total electron scat- 
tering optical depth r = O.OSS^jJ I]^^ (Spergel et al. 2005, Table 
5) that has been inferred from the WMAP three-year polariza- 
tion data (Page et al. 2005) implies that reionization occurred 
relatively recently, at Zr ^ 1 1 ■ ' 

An outstanding question, and the focus of this Letter, is whether 
the gas in the majority of the first-generation low-mass halos 
was able to cool and form ionizing sources (stars and/or black 
holes [BHs]), and contribute to reionizing the IGM. The neu- 
tral gas in halos with virial temperatures Tvii ^ lO'^K (hereafter 
"minihalos") cannot cool via atomic hydrogen alone, and the 
presence of H2 (or other) molecules is a necessary condition 
for the gas to reach high densities and to ultimately form stars 
or BHs in these halos (Saslaw & Zipoy 1967). As a result, these 
halos are distinct from larger halos with J> 10'*K, which 
are expected to have different star-formation properties (Oh & 
Haiman 2002; Johnson & Bromm 2006). Since the smallest ha- 
los are the first to collapse, whether or not typical minihalos are 
able to host ionizing sources is the most significant factor in de- 
termining when reionization began, and driving the value of r 
(Haiman & Holder 2003 [HH03]; see Abel & Haiman 2000 for 
a review focusing on the role of H2 molecules for reionization). 

In isolation, minihalos with virial temperatures as low as a 
few 100 K could form enough H2, via gas-phase chemistry, for 
efficient cooling and gas contraction (Haiman, Thoul & Loeb 
1996; Tegmarket al. 1997). However, H2 molecules are fragile, 
and can be dissociated by soft UV radiation absorbed in their 

' Assuming an abrupt transition from a fully neutral to a fully ionized IGM, 
and assuming helium is once ionized at the same redshift as hydrogen, we find 

Zr < 10. 

(2005). 



Lyman-Werner (LW) bands (e.g. Haiman, Rees & Loeb 1997; 
Ciardi et al 2000; Ricotti, Gnedin & ShuU 2001). In patches 
of the IGM corresponding to fossil HII regions that have re- 
combined, the gas retains excess entropy, which can reduce gas 
densities in the cores of collapsing halos, and lower the LW 
radiation background that photodissociates H2 (Oh & Haiman 
2002; Gnedin 2000). On the other hand, positive feedback ef- 
fects, such as the presence of extra free electrons (beyond the 
residual electrons from the recombination epoch) from proto- 
galactic shocks (Shapiro & Kang 1987; Ferrara 1998), from 
a previous ionization epoch (Oh & Haiman 2003; Susa et al. 
1998), or from X-rays (Haiman, Rees & Loeb 1996, Oh 2001; 
Ricotti et al. 2002a,b), can enhance the H2 abundance. Whether 
or not gas cooling in minihalos was efficiently quenched glob- 
ally has remained unclear, with numerical simulations generally 
favoring less quenching (Machacek et al. 2001, 2003; Ricotti 
et al. 2002a,b; Kuhlen & Madau 2005) than predicted in semi- 
analytical models. 

The purpose of this Letter is to place constraints on the ef- 
ficiency of minihalos to produce and inject ionizing photons 
into the IGM, using the value of r inferred from the WMAP 
polarization data. Although the absolute values of the relevant 
efficiency parameters are uncertain, there is a robust lower limit 
on their combination from the requirement that reionization is 
completed prior to z = 6. As a result, a measurement of r per- 
mits robust conclusions about the evolution of the efficiency 
with redshift or halo mass-scale (HH03; Cen 2003; Onken & 
Miralda-Escude 2004). The main conclusion of this Letter is 
that the efficiency with which minihalos injected ionizing ra- 
diation into the IGM was about an order of magnitude lower 
than expected and lower than in larger halos, requiring negative 
feedback at high redshift. 

Throughout this letter, we adopt the background cosmologi- 
cal parameters as measured recently by the WMAP experiment 
(Spergel et al. 2005, Table 5), Vl,„ = 0.24, VLa = 0.76, Vlb = 
0.0407, /z = 0.72 and an initial matter power spectrum P{k) oc 
with n = 0.96 and normalization (Tg = 0.76. The values of these 
parameters determine the required overall efficiencies, but do 
not otherwise have a significant effect on our conclusions. 
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2. MODELS OF REIONIZATION 

In this section, we briefly describe our semi-analytical model 
of the reionization process. The treatment follows the prescrip- 
tions in HH03, and the reader is referred to this paper for a 
detailed description; we only recapitulate the main features and 
a few important modifications here. The models track the total 
volume-filling factor Fuu(z) of ionized regions, assuming that 
ionizing sources are located inside virialized dark matter ha- 
los. The sources create ionized HII regions, which expand into 
the IGM at a rate dictated by the source luminosity, and by the 
background gas density and clumping factor^ 

The efficiency with which the ionizing source(s) associated 
with each halo inject ionizing photons into the IGM is parame- 
terized by the product e=Nyf^,f esc, where/* =M*/(ilbMhaio/f^m) 
is the fraction of baryons in the halo that turns into stars; A^-^ 
is the mean number of ionizing photons produced by an atom 
cycled through stars, averaged over the initial mass function 
(IMF) of the stars; and /esc is the fraction of these ionizing 
photons that escapes into the IGM. We define separate efficien- 
cies, Emini, fiaige for the minihalos and large halos, respectively. 
The distinction between the halos is made at Tyir = lO'^K, with 
the relation 7;i,- « 18OO(M/1O''M0)2/3(1 +z)/21 K. We adopt 
the virial temperature for the smallest minihalos that can cool 
via H2 and form stars to be Tj^in = 400K (Haiman et al. 1996; 
Tegmark et al. 2001; Machacek et al. 2001 - note that the 
choice in HH03 was lower, r^in = 100 K). 

Gas infall is expected to be suppressed for sufficiently small 
halos in actively ionized and photo-heated patches of the IGM. 
We define a temperature and only permit halos with virial 
temperatures below this value to form ionizing sources in the 
neutral fraction 1 - fHii(z) of the volume of the IGM. The value 
of Tuv is uncertain - at low redshifts it is commonly taken as « 
2.5 X 10^ K, corresponding to a velocity dispersion cr = 50 km s"' 
(e.g., Thoul & Weinberg 1996). However, Dijkstra et al. (2004) 
showed that this suppression is less efficient at high redshift, 
and based on their result for halos collapsing at redshifts 6 ^ z ^ 1 0, 
we adopt the fiducial value of Tuv = 9 x lO'^K (or a ^ 30km s"'). 

In order to match inferences from the z ^ 6.5 quasar spectra 
(Becker et al. 2001; Fan et al. 2006), we require that percola- 
tion takes place at redshift z = 6.5. At this redshift, ionization is 
dominated by larger halos, and so this requirement effectively 
fixes the value of eiarse = 120. This value of the overall ef- 
ficiency is reasonable: for a normal stellar population with a 
Salpeter mass function, Nj = 4000, and we may break down 
the efficiency into /, = 0.15 and /esc = 0.2. The star-formation 
rate at z ^ 6 predicted with these assumptions is consistent with 
the global value inferred from the Hubble Ultra Deep Field (see 
Mesinger et al. 2005 for a comparison). 

Finally, the evolution of the ionized fraction depends on the 
clumping factor of the ionized gas, Chii = («Hii)/("Hn)^- The 
value of the clumping factor at z 6 does not impact our re- 
sults because it is absorbed into eiarge, which, in turn, is set by 
the requirement that the ionization percolates at z = 6.5. This 
normalization requires eiarge ~ 120 x Chu(z = 6)/10. However, 
we find that the evolution of Chii with redshift has a significant 

^ We do not explicitly treat the reionization of helium, but assume the Hell/Hel 
fraction follows the HII/HI fraction when we compute r. For neutral/doubly 
ionized helium, r would be smaller/larger by ~ 8%. 

^ HH03 found a value that was 50% smaller - we here find that the efficiency 
increases by a factor of 3 because of the lower small scale power impHed by the 
3-yr WMAP data (lower erg and «.,), but was decreased by a factor of 2 due to 
our lower value of Tuy . 
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Fig. 1 . — The assumed evolution of the clumping factor of ionized gas with 
redshift. The thick solid (red) curve is adopted as our fiducial case, and has 
a decrease by a factor of ^ 4 between z = 6 and z = 15. The shaded range 
indicates models in which this decrease varies from a factor of 2— 10; the 
hmiting case of a constant clumping factor of Chii = 10 is shown as the dotted 
hne. 



effect on our results. We do not attempt to model this evolu- 
tion ab-initio; instead we adopt a relation that is a conservative 
compromise between three different ab-initio approaches, and 
we vary our parameters to reflect the range of predictions of the 
three models (see discussion in § 3 below). The relation we 
adopt is a simple power-law, 

Cnn(z)-1 = 9(^^^Y , (1) 

which has a fixed Chii(z = 6) = 10, and declines towards high 
redshift at a rate given by the slope (3. What matters most for 
our results below is the decrease between z = 6 (the redshift 
which sets eiarge) and 15 z ^ 20 (the redshift which determines 
the contribution to t from the minihalos). In our fiducial model, 
13 = 2, coiTesponding to the ratio Chii(z = 6)/Chii(z = 15) = 3.7. 
This choice, together with the range 1 < /3 < 6 (corresponding 
to a range of ratios 2 < Chii(z = 6)/Chii(z= 15) < 10), is shown 
in Figure 1 . For reference, the dotted line shows the limiting 
case of a constant clumping factor (/3 = 0). 

In summary, our models have only two free parameters: emini, 
and (3. In addition, we will consider uncertainties in Tmin, the 
minimum mass for efficient H2 cooling, and Tuv, the maximum 
virial temperature of halos that can collect photo-heated gas. 
We next discuss our choices for the fiducial values of these pa- 
rameters. 

3. FIDUCIAL PARAIVIETERS 

What is the expected value of emmi? Abel, Bryan, & Norman 
(2002), and Bromm, Coppi & Larson (2002) suggest that mini- 
halos form metal-free stars from only a small fraction (/* « 
0.0025) of the available gas. On the other hand, the stars are, on 
average, w 20 times more efficient ionizing photon producers 
(per unit stellar mass) because they are massive (M ^ 100 Mq) 
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Fig. 2. — The evolution of the ionized fraction of hydrogen. The thick .solid 
(red) curves, the shaded ranges, and the limiting cases shown by the dotted 
curves, correspond to the different logarithmic slopes for the evolution of the 
clumping factor with redshift Chii(z) shown in Figure 1. The lower set of four 
curves assume star-formation in minihalos (with 400K < Tyi^- < Iff^K) has 
been completely suppressed (£,„!„] = 0), and the upper set of curves assumes 
minihalos inject ionizing photons into the IGM with the fiducial efficiency 
(e™„i = 200). 



and have harder spectra (Tumhnson & Shull 2000; Bromm, Ku- 
dritzki & Loeb 2001; Schaerer 2002). As a result, A^..^ « 80000; 
and all of thek ionizing radiation escapes into the IGM, at least 
from the smaller minihalos that contribute most of the high- 
redshift tail for reionization (/esc = 1; Whalen, Abel & Nor- 
man 2004). This results in the overall fiducial efficiency of 
emini = Cfid = 200. We note that this is quite close to the value 
eiarge =120 that we inferred for the large halos by requiring per- 
colation to occur at z w 6.5. 

Several authors have attempted to compute the clumping fac- 
tor Chii(z) from first principles. In the simulations of Gnedin & 
Ostriker (1997), it increases dramatically (by a factor of ^ 50) 
between z = 6 and z = 15, but since they do not resolve small- 
scale fluctuations early on, this ratio may be an overestimate. 
Iliev, Scannapieco & Shapiro (2005) have a different approach 
of separately estimating the contribution from the low-density 
IGM from a simulation, and adding to this semi-analytic es- 
timates of the clumping due to minihalos (Haiman, Abel & 
Madau 2001). They find an evolution that is less steep, with 
Chii(z =10) /Chii(z = 15) w 2, but they do not cover the interest- 
ing z < 10 range, where we normalize our models, and where 
Chii likely increases rapidly. Finally, Miralda-Escude, Haehnelt 
& Rees (2000) offer a different prescription: they assume that 
the IGM is ionized "outside-in", with all the gas below a crit- 
ical overdensity Acrit ionized at a given redshift (and neutral 
above this threshold), to compute Chii directly. This prescrip- 
tion predicts relatively slow evolution, with Acrit ~ 10-20, un- 
til the late stages of reionization (as in Wyithe & Loeb 2003 
and Wyithe & Cen 2006). However, towards the end of the 
reionization epoch, Acrit increases rapidly. This is an important 
distinction for our present purposes, because the ionizing emis- 
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FlG. 3. — The optical depth to electron scattering t, as a function of the 
efficiency of minihalos to inject ionizing photons into the IGM. The efficiency 
is shown in units of the fiducial value ef^ = 200 (see text for discussion). The 
thick solid (red) curve, the shaded range, and the limiting case shown by the 
dotted curve, follow the notation for models with different clumping factor 
evolutions shown in Figures 1 and 2. The horizontal lanes show the range of t 
allowed by the WMAP three-year data (at 68% and 95% CL). 



sivity in our models is normalized at the end-stages of reion- 
ization. Miralda-Escude, Haehnelt & Rees (2000) show that by 
the time the mean free path of ionizing photons grows to the 
Hubble length (which appears to be the case at z ^ 6; Lidz et 
al. 2005), Acrit 500, which will correspond to a steep rise in 
the clumping factor. 

4. RESULTS AND DISCUSSION 

The evolution of the ionized fraction of hydrogen in several 
models is shown in Figure 2. The thick solid (red) curves, 
the shaded ranges, and the limiting cases shown by the dotted 
curves, correspond to the different clumping factor evolutions 
Chii(z) shown in Figure 1 with the same notation. The lower set 
of curves assumes star-formation in minihalos has been com- 
pletely suppressed (emini = 0), and the upper set of curves as- 
sumes minihalos inject ionizing photons into the IGM with the 
fiducial efficiency (emini = 200). All curves are normalized to 
produce full reionization at z ^ 6.5. Earlier reionization is con- 
sistent with SDSS quasar spectra; this would yield a larger r, 
and would strengthen our conclusions. 

The models shown in Figure 2 that include minihalos have 
optical depths of r = 0.14,0.17,0.19, and 0.23, and are all ruled 
out by the three-yeai" WMAP data at more than 95% CL. The 
models that exclude minihalos have t = 0.077, 0.084, 0.091, 
and 0.12, and all except the case with the steepest Chii(z) evo- 
lution are within 1 cr of the WMAP value. 

In Figure 3, we show t as a function of the efficiency Cmini. 
The figure follows the notation of Figures 1 and 2 to describe 
the four different cases for the evolution of the clumping factor. 
The main conclusion that can be drawn from this figure is that 
the photon production efficiency in minihalos was smaller than 
the expected fiducial value by a factor of ^2-10. This also 
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Fig. 4. — The predicted optical depth as a function of the ionizing efficiency 
of the minihalos (emini/ffld) and the logarithmic slope of the decline of the 
clumping factor with redshift (/3 = fidogCCmi - l)/(/log(l +z)). The lowest 
contour (white region) corresponds to optical depths within one sigma of the 
3-yr WMAP result, and each contour is a one-sigma increase in the optical 
depth. The fiducial model is marked with an asterisk and deviates from the 
WMAP value by nearly 4a. 



means that minihalos contribute to reionization with a lower 
efficiency than the larger halos. 

More generally, Figure 4 shows the predicted optical depth 
as a function of both €,mm and (3. The lowest contour (white 
region) corresponds to optical depths within Icr of the 3-yr 
WMAP value, and each successive contour shows a 1<t increase. 
The figure shows that our fiducial model deviates from the WMAP 
value by nearly 4a; consistency at the I - 2a level requires ei- 
ther reducing emini by a factor of 2-10, or adopting an essen- 
tially non-evolving clumping factor (/? « 0) and a smaller re- 
duction in Cmini- 

In order to examine the robustness of our conclusions to un- 
certainties in the model, we have repeated our calculations for 
different values of T„^n, the virial temperature of the smallest 
minihalo, and Tuv, the critical virial temperature for suppression 
of gas infall by photoionization heating. We find that variations 
in the range 160 K < r,ni„ < 10^ K and 10"^ K < T^v < 2 x 10^ K 
change the required reduction factor of Cmini by less than 30%. 

We emphasize that our predicted range of r = 0.08-0.11, 
for models that exclude minihalos, is consistent with numerous 
previous semi-analytic studies (as well as simulations that do 
not resolve minihalos, e.g., Ciardi et al. 2003), in the wake 
of the 1-yr announcement of the WMAP measurements. In 
another recent study, Wyithe & Cen (2006) examined popu- 
lations of metal-free stars in Tyir <; 2 x 10'*K halos, with an 
efficiency similar to normal stars, and found similar r values. 
It is also worth noting that the modified Press-Schechter form 
for the DM halo mass function that we used (Jenkins et al. 
2001) has been directly probed and confirmed in the relevant 
mass and redshift range by recent three-dimensional simula- 
tions (Yoshida et al. 2003; Jang-Condell & Hernquist 2001; 
Springel et al. 2005). However, the minihalos arise from pri- 
mordial perturbations on very small (A: ~ 10- 100) scales, and 
predicting their abundance does involve an extrapolation of the 
measured power spectrum by -^2-3 orders of magnitude in k 
(Barkana et al. 2001; Mesinger et al. 2005). 



5. CONCLUSIONS 

The WMAP experiment has opened a new window into stud- 
ies of the first structures at the end of the cosmological dark 
ages. The first-year data release has suggested the large optical 
depth T ~ 0. 17 (Spergel et al. 2003), which appeared in conflict 
with simple predictions, and resulted in different suggestions to 
account for this large r. We argue here that the significantly 
lower r « 0.09 ± 0.03 allowed by the three-year data is once 
again interesting: if minihalos were forming stars efficiently, 
the expected value would be higher, r ^ 0.17. To reproduce the 
new value for t requires an order of magnitude reduction in the 
fiducial efficiency of ionizing photon production in minihalos. 

Star-formation in minihalos has been predicted theoretically 
to be significantly suppressed, and the new WMAP data has pro- 
vided the first interesting empirical evidence for this suppres- 
sion. In addition to the various physical feedback effects men- 
tioned in the Introduction, the clustering of the early minihalos 
will, in general, help suppressing more halos at high redshift, 
since clustering will place more minihalos in the destructive 
sphere of influence of earlier collapsed objects. In the presence 
of clustering, photoionization heating iilone may account for 
the low r (Kramer et al. 2006). 

If massive Pop HI stars typically leave behind intermediate- 
mass black holes, these may accrete efficiently and produce 
more significant ionization than their progenitor stars (Madau 
et al. 2004; Ricotti & Ostriker 2004). Our results, by inference, 
provides similarly interesting constraints on this scenario. 

Further studies from CMB anisotropics, spectroscopic ob- 
servations of quasars and Lya emitting galaxies at z > 6, radio 
probes of the redshifted 21cm fine of neutral hydrogen, as well 
as direct detections of the sources of reionization and their end- 
products as supemovae and gamma ray bursts, wiU clarify the 
star-formation efficiency in minihalos and their contribution to 
reionization. 
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